On the low fermionic eigenmode dominance 
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We demonstrate the utility of a spectral approximation 
to fermion loop operators using low-lying eigenmodes of the 
hermitian Dirac- Wilson matrix, Q = 75M. The investigation 
is based on a total of 400 full QCD vacuum configurations, 
with two degenerate flavors of dynamical Wilson fermions at 
(3 = 5.6, at two different sea quark masses. The spectral 
approach is highly competitive for accessing both topologi- 
cal charge and disconnected diagrams, on large lattices and 
small quark masses. We propose suitable partial summation 
techniques that provide sufficient saturation for estimating 
Tr , which is related to the topological charge. In the 
effective mass plot of the ri' meson we achieved a consistent 
early plateau formation, by ground state projecting the con- 
nected piece of its propagator. 

PACS: ll.15.Ha, 12.38.Gc 



I. INTRODUCTION 

The Euclidean approach of lattice gauge theory (LGT) 
has been established as a viable framework to deal with 
quantum chromodynamics beyond the limitations of per- 
turbation theory [Q. Whereas practical algorithms have 
been developed to calculate many physical observables, 
there is a class of physical quantities that is highly sensi- 
tive to fluctuations in the vacuum gauge field that has not 
been amenable to known ab initio numerical methods. 
Notorious examples of such computationally intractable 
problems include mixing phenomena between fermionic 
and glueball states dj, disconnected quark loop diagrams 
occurring in flavor singlet matrix elements , and the 
infamous rj' propagator [||-|7|. 

In all these instances one would like to take advantage 
of self-averaging effects by exploiting the translational 
invariance of the QCD ground state. This amounts to 
probing with objects shifted across all lattice space-time 
points and hence to the computation of light quark prop- 
agators on the entire vacuum field, i.e. the full inverse of 
the Dirac operator, M, which exhibits a high condition 
number in the regime of light quark masses. Unfortu- 
nately this is a prohibitively expensive numerical task, 
with M being of rank > 10^ on typical lattice sizes. 



Therefore the standard approach is to resort to stochastic 
estimator techniques [p[-p^ . 

The approach to these previously intractable prob- 
lems explored in this work is motivated by the expec- 
tation that low-lying modes of the Dirac operator should 
embody the important features of fermionic physics in 
the chiral regime Q Physically, topological excita- 
tions in the QCD vacuum, corresponding to instantons in 
the semiclassical limit, generate low-lying fermion modes 
that play an essential role in light quark physics, rang- 
ing from generating the 't Hooft interaction to producing 
the chiral condensate. Hence, for sufficiently light quark 
masses, expansion in a basis that contains these low-lying 
eigenmodes should cover the essential physics associated 
with observables involving the topological charge, fiavor 
singlet disconnected diagrams, and processes described 
by the 't Hooft interaction. 

It therefore appears worthwhile to launch another 
attempt to explore the potential of spectral meth- 
ods/approximations [ p^ . The question is, how many 
low-lying eigenmodes suffice to bear out the important 
features of long range physics in practical instances. 

There are two options to proceed, based on the spectral 
representations of the non-normal Wilson-Dirac matrix 
M and the hermitian matrix Q — "f^M . The advantage of 
M lies in its shiftable (with respect to k) structure M = 
1 — kD, but it requires to work with biorthogonal sets 
of eigenmodes and complex eigenvalues. The hermiticity 
of Q, on the other hand, allows for a simple ordering 
of the orthogonal eigenmodes and a natural definition of 
low-lying eigenmodes. 

We decided to work with Q := 75M and its eigen- 
modes [EMll 
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We can express the quark propagator in terms of the 
'eigenmodes' (Ai, iV'i)) from the spectral representation: 



Q ^(n, a, a; m, b) — 
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^ Early pionee ring w ork in this direction has been done by 
Barbour et al. @l|l. 
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FIG. 1. Distribution of moduli of the eigenvalues of Q, av- 
eraged over the configurations. The upper line corresponds 
to K = 0.156, the lower to k = 0.1575. 
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FIG. 2. Same as Fig. but plotted in terms of eigen- 
value ratios, |Ai(0.1575)/Ai(0.1560)|, both ordered according 
to Eq. 0. 



where n, a, a etc. label lattice sites, Dirac and color in- 
dices, respectively. 

For our practical benchmarking we use 200 (195) 
SESAM lattices of size 16^ x 32 with k = 0.1575 
(0.1560) ||. Thus the rank of Q is dim(Q) = 1 572 864, 
which bars us from complete diagonalization of Q and 
the use of the identity as such, Eq. (||). Nevertheless, the 
low-lying eigenvalues, i.e. the ones with smallest moduli, 
represent a large weight and are expected to saturate the 
sum in the regime of small quark masses. 

In the present paper we wish to investigate to what 
extent the spectral sum allows, under the conditions of 
the SESAM simulations, for a truncation of the propaga- 
tor at a reasonable (in the sense of practicality) number 
of lowest-lying eigenmodes. As a testbed for our trun- 
cated eigenmode approach (TEA), we will use Tr in 
addition to the pion and r/'-correlators. 
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FIG. 3. Partial sums, according to Eqs. and from a 
particular SESAM configuration at k = 0.1575. 



II. COMPUTING LOW-LYING EIGENMODES 

The hermitian operator Q is defined through 

Q = 15M , 

where M is the standard Wilson-Dirac matrix 

M{n, a, o; m, (3, b) — ——S{n, a, a; m, (3, b) (3) 
2k 

^ 4 

~2 ^ [('^(";/5) - Ip-ia; P))Us^^f,{a;b)d{n + fi;m) 

+ {S{a-p)+j^{a;p))ul_^Ja;b)S{n~fi; m) 

Here the 7^ are the Dirac matrices, U is the gauge field 
(vacuum configuration), d is the Kronecker delta function 
and p, denotes a unit vector in a space or time direction. 

The eigenmodes are determined by use of the paral- 
lelized Arnoldi package^ . We speed up the computa- 
tions by application of the polynomial acceleration tech- 
nique [ P4[ . 

The Arnoldi method works efficiently when calculating 
eigenvalues on the surface of the spectrum, here around 
Amm and Xmax- For this reason we first need to execute 
a preparatory step where we map, by means of a suitable 
polynomial q : Q t— > q{Q), the low-lying modes from the 
inside of the spectrum to the surface. We will contrive the 
polynomial such that the spectrum is already prepared 
for the subsequent Chebyshev acceleration step |2^ . 



^We employ the Arnoldi method as provided by the parallel 
Arnoldi Package PARPACK Q from Rice University. The 
Arnoldi method is designed for non-hermitian matrices, but it 
reduces, when applied to a hermitian matrix, to the Lanczos 
method. 
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FIG. 4. Same as in Fig. |^. The full lines connect the entries, t-' , with equal values of j. The horizontal line corresponds to 
the value t as defined in Eq. ^ 



a. First step: Mapping We apply a very simple poly- 
nomial transformation q. To be specific we proceed as 
follows: 

- Target a spectral window of Q within which all 
eigenmodes are to be determined and decompose 
the spectrum of Q accordingly into desired eigen- 
modes, 



W = spec(Q) 



(4) 



and the remainder: 

spec((3) = spec{Q)w U spec((3)„ . (5) 

Construct a polynomial q that casts the spectral 
window spec{Q)w outside the region [—1, 1]: 

{q{X)\X e spec(Q)„} C [-0O, -1[ U ]1, oo] , (6) 
{(?(A)|Aespec(Q)4c [-1,1] . (7) 



We found that the simple choice 
2 



qiQ) 



;Q'-il + r)I 



(8) 



will accomplish the task, since the smallest eigenmodes 
ofq{Q), i.e. the ones closest to Xmin{q{Q)), correspond to 
the low-lying eigenmodes of Q. The polynomial carries 
two parameters: 

(i) The scale factor s is the spectral radius of Q- It can 
be computed in a first Arnold! run on just a few vacuum 
configurations, since it fluctuates little with the gauge 
field {[/}. 

(ii) The offset parameter r represents a simple shift oper- 
ation and controls the actual size of the window of desired 
eigenmodes. 

b. Second step: polynomial acceleration The Arnoldi 
method is sensitive to the level density of eigenvalues, 
p(A). For acceleration it is therefore crucial to precondi- 
tion the problem by decreasing the level density in the 
spectral region of interest. This can readily be achieved 
with Chebyshev polynomials of degree TV due to their 
rapid increase outside the interval [-1,1] (within which 
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FIG. 5. Partial sum, showing for each s{l) the associated j 
such that s{l) = t^l). 



they are close to zero). 

The practical procedure is then to find an appropriate 
Chebyshev polynomial and compute the eigenmodes 
of largest moduli of Tn o q o Q, with q taken from Eq. (H). 

We found the optimal convergence w.r.t. CPU time 
with a Chebyshev polynomial Tn of order = 80 and 
r chosen such that approximately only the searched for 
300 lowest-lying eigenvalues were contained in the inter- 
val [— 1 — r, — 1]. With this parameter setting we needed 
just one Arnoldi factorization Since we chose the 

size of the factorization to be 600, the calculation of the 
eigenmodes takes 2 x 80 x 600 = 96000 matrix-vector 
multiplications (where the factor 2 is due to the fact that 
enters the Chebyshev polynomial). 

We emphasize that without the Chebyshev accelera- 
tion step, we observed only very poor convergence. 

Once we have the eigenf unctions, \ipi), of the operator 
Tn ° q ° Q, it is straightforward to retrieve the eigenval- 
ues Xi of the original operator Q. By appealing to the 
relation 

QIV-*) = A,|V^) ^ TN{q{Q)m) = TN{q{X^m^) , (9) 

we determine through the Rayleigh quotient 



(10) 



Let us comment that - for our present purposes - we 
found the above procedure to be much superior to the 
shift and invert strategy which is frequently used to 
find the low-lying eigenvalues 

In Fig. ([l]) we illustrate the spectrum of Q averaged 
over the SESAM configurations, at lightest and heaviest 
sea quark mass respectively, for which we plot the de- 
pendence of |Ai| on i. The ordering is chosen according 
to 
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FIG. 6. High resolution plot of with 600 Eigenvalues 



|Ai| < IA2I < IA3I < 



(11) 



The upper (lower) line corresponds to the heaviest (light- 
est) sea quark mass of the SESAM sample. Their ratio, 
I Ai(hghtest)/Ai (heaviest) I is plotted in Fig. (^) as a func- 
tion of i. It shows a rather steep rise for i < 100, which 
illustrates the growing importance of the lowest-lying 
modes in the spectral representation of Q~^, Eq. (||), 
when decreasing the sea quark mass. This feature should 
become even more pronounced for the imminent QCD 
simulations in the yet deeper chiral regime, mps/mv < 
0.5. 



III. DETERMINATION OF Tr Q"^ 

In the continuum limit, the complete trace 

S = Tr (12) 

is related via the index theorem |^ to the net topo- 
logical charge of the gauge field. Albeit the breaking of 
chiral symmetry by the standard Wilson lattice action 
blurs this connection somewhat at finite lattice spacings, 
a, one would still expect remnants of the index theorem 
to hold ll^. In previous works, Q has been estimated by 
use of stochastic estimator techniques (SET) which basi- 
cally compute quark propagators on stochastic sources, 
and indeed, a close correlation between gauge field topol- 
ogy and Tr Q"^ has been observed |p9| , pO| . Hence, it is a 
quantity of significant interest. This will be all the more 
the case in future applications of the overlap fermion lat- 
tice formulation |3l|| . 

Q is given by the sum of the inverse eigenvalues: 



TrQ- 



E 



1 

Xi {'^l^i\'lp^) ^ A,; 



(13) 
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FIG. 7. Functions as in Fig. n for 10 gauge fields such that their t-values are all equal to 0. 



Physically, it is the low-lying modes that encapsu- 
late the interesting information within the full sum of 
Eq. (|l^) while the large modes provide just nuisance by 
adding background noise to the infrared signal. Heu me 
miserum, in order to turn this qualitative proposition 
into a quantitative statement we would need to know the 
transition point between infrared and ultraviolet physics. 
For an approximate use of the spectral relation, Eq. (13), 
we need a cunning technique to deplete the unwieldy 
background. 

With this in mind let us take a heuristic approach and 
study the pattern of eigenvalue distributions as obtained 
from the SESAM ensemble of vacuum configurations. To 
reach our goal, it is useful to order the spectrum and 
perform certain partial summations in Eq. (p^): 

Let Pi denote the positive and rii the negative eigen- 
values, 



{AJ = {p,} U {nj, 
and let them be ordered such that: 



Pi <P2 < P3 



< . . 



\ni\ < \n2\ < Inal < . . 



(14) 

(15) 
(16) 



In order to study convergence properties let us define 
the following partial series: 



j+k 



j+k k 



J >o 

j < 



l^\j\ + 2k, 



(17) 

(18) 
(19) 



where, in obvious notation, the parameter j labels the 
excess number of entries with positive (negative) over the 
ones with negative (positive) eigenvalues. In Fig. (|^) and 
Fig. (H) the values of P are plotted for the index range 
~5 < J < 5, as obtained from one particular SESAM 
configuration at our lightest available quark mass, k = 
0.1575. The partial sums appear to exhibit a certain 
convergence pattern which is displayed in Fig. (^, where 
we connected the points to given j-values. 

Alternatively, one might organize the summation in the 
order of increasing moduli of the eigenvalues (Eq. (pi])), 
independent of their signs, and define the truncated sums 
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FIG. 8. Functions as in Fig. ^, but for a heavier quark mass. 
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(20) 



The family of curves, J- — {^■'(0} provides a suitable 
framework to disclose the asymptotic behavior of this 
inverse eigenvalue summation, s{l). This is illustrated, 
again for the particular gauge configuration, in Fig. (j^), 
where we display the data points for s{l) in the range 
1 < Z < 300. For reference we also show the (slowly 
narrowing) band of the partial sums, Note that 

from s ~ 150 onward, s{l) in this particular configuration 
jump^mostly between the two levels and t~'^. More- 
over t appears to be distinguished as lying between the 
asymptotically falling set of curves with j > — 2 and the 
rising ones, j < —4. 

We found this scenario to apply mutatis mutandis to 



^ To put it differently: the asymptotic distribution of eigen- 
values is characterized by alternating signs when proceeding 
according to the order given by Eq. (|ll|). 



all configurations in the sense that 

• for each gauge field {U}, T contracts around a par- 
ticular partial sum tP{l) that levels to a plateau 
value beyond I « 150, with p depending on the 
choice of the gauge field [U]. 

Let us next quantify our observations and denote the 
height of this supposed plateau with t, setting 



t := iP(300) 



(21) 



Note that t varies with the underlying gauge configura- 
tion [U] . The question then arises how accurately we can 
extract the actual plateau height. In order to appreciate 
the numerical flatness of the plateau curve we pushed the 
eigenvalue computation to a number of 600 modes on a 
single SESAM configuration. In Fig. ^ we display the 
resulting plateau on a magnified scale. Assuming that 
the apparent remaining weak oscillation for I > 150 is a 
pars pro toto feature for the entire SESAM sample and 
that with 300 eigenvalues one has already passed the first 
extremum on the entire sample of gauge fields J7, we es- 
timate that from the lowest-lying 300 eigenvalues t can 
be determined with a bona fide accuracy of 1%. 
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FIG. 9. Pattern of partial series t-' on a quenched 4* lat- 
tice at P — 5.0, which allows for a calculation of the entire 
spectrum. 
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Tr Q"' from TEA 

FIG. 10. Scatterplot comparing the results from TEA and 
SET on all configurations at k = 0.1575. The errorband is 
chosen such that 67% of the data points lie within. 



Next let us argue that the plateau value, t[U], provides 
us with an approximant for the complete trace, Eq. (12): 



Q[U] w t[U] 



(22) 



What is the deviation from the complete trace in our sit- 
uation? In our range of quark masses, we can exclude 
zero-level crossings of eigenmodes. Therefore the matrix 
Q possesses an equal number of positive and negative 
eigenvalues. Hence, when adding up all n eigenvalues 
(with n — dim((5)), s{n) will lie on the curve with super- 
script 



Q = s{n) = (n) , with n = dim(Q) . 



(23) 



Thus Tr Q ^ is related to the plateau height t in the 
following way: 



Q = t 



P 



(24) 



FIG. 11. Comparison of tea's estimates of the topological 
charge with the gluonic determination after cooling, on all 
SESAM configurations at lightest quark mass k — 0.1575. 



The second term on the r.h.s. measures the distance be- 
tween and tP at I = n. It can be neglected with respect 
to the error on t itself, since 1 1 / A„ « . 1 5 1 is approximately 
equal to the error of t. 

Another justification for the validity of our approxima- 
tion, Eq. ( p^ , comes from our observation that different 
field configurations [Ur] yield equal results when plotted 
with appropriate offsets t[Ur], namely (P {l)[Ur] ~t[Ur])- 
This is illustrated in Fig. (0) and Fig. (||) where we su- 
perimpose 10 such series for Ur{r = 1, . . . 10): it strikes 
the eye that, for I > 150, the partial sums to the 10 
gauge fields all collapse onto a single, universal family of 
curves: 



t^'+f -(/)[[/,] - t[Ur] = P+P^^ il)[U,] ~ t[U, 



(25) 



Hence there exists a set of [/-independent functions t^{l) 
such that the following identity with respect to U applies: 



ti+PVJ\(l)[u]-t[U]=P{l) 



(26) 



This pattern strongly supports the picture that on our 
configurations the interesting physics with respect to the 
topological charge is indeed contained in the subset of 
the 150 smallest eigenvalues, while the remaining ones 
carry no information on Q. 

We corroborate this result by considering a 4* lat- 
tice where we determined all 3072 eigenvalues of Q in 
quenched QCD at /3 = 5.0. In Fig. (|^) we plotted the 
corresponding partial sums P . Notice that V (here p 
happens to be 0) remains absolutely fiat after reaching 
its plateau value at around I ~ 500. 

A comparison of the Q values as produced by TEA 
with the ones obtained in Ref. |3^ from SET on the en- 
tire K = 0.1575 sample is shown in Fig. (|l0|). The data 
points seem to scatter rather nicely around the bisecting 
line. While the accuracy of the TEA results on individual 
configurations is about 1%, the uncertainty of the SET 
estimates turns out to be AQ « ±50. 
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FIG. 12. Set of TT correlation function C'{At) from TEA 
on local sinks and sources, plotted against the spectral cutoff, 
I. At increases as one steps down from the top curve which 
refers to At — 0. 
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FIG. 13. Comparison of the n correlation function as pro- 
vided by TEA with / — 300, with the one obtained from the 
standard method (solving linear systems) on local sinks and 
sources. 



In Fig. (grj) we present TEA's (normalized) Q values 
along the Monte Carlo history of our SESAM sample 
at K = 0.1575 and compare them to the result of the 
gluonic determination after cooling |^^. We reconfirm 
our previous finding ]29| ] that there is a close correlation 
between the gluonic and fermionic definitions of the net 
topological charge. 



FIG. 14. Local tt masses from TEA and the standard 
method . 
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FIG. 15. Local masses of the tt from TEA. LL are the un- 
smeared results, whereas SS stands for smearing the sinks and 
sources. The horizontal lines show the errorbands of the tt as 
obtained from the standard method (solving linear systems). 



ular hadronic channels. 

In order to set the stage we shall first consider the 
octet pseudoscalar channel as we can easily compare to 
standard pion correlator computations. We shall then 
elaborate on the spectral approach to the singlet pseu- 
doscalar propagator which differs from the octet one by 
the two-loop correlator. The latter involves the compu- 
tation of incomplete trace expressions {Q{t)Q{t + At)) 
and thus represents a quantity of increased complexity. 



IV. HADRONIC TWO-POINT FUNCTIONS 

In this section we wish to investigate the potential 
of spectral methods in the computation of two-point 
hadronic correlators j32). The question here is to what 
extent we can verify low-lying eigenmode dominance for 
infrared physics, specifically the ground states in partic- 



A. Basics for mass determinations 

Masses are extracted from the large time behavior of 
correlation functions C{t) that carry the quantum num- 
bers of the particles in question. The correlation func- 
tions of the TT and the rj' have the following form 
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FIG. 16. The dependence of the local masses m{At) on 
the spectral cutoff, I < Imax ~ 600, on a single configura- 
tion. The X-axis carries the parameter l^t ~ 600At -f / with 
At — 1,2,...,15. The horizontal line shows the errorband of 
the TT mass as obtained from the standard method (solving 
linear systems). 

C^(At) = ^ ^ [Q~'^{si,t,ai,ai;s2,t + At,a2,a2) 
Q~^{s2,t + At, a2,a2] si,t, ai, ai)] (27) 

c^>{At) = c^{At) - ^/( I] Q{t)Q{t + 

with Q{t) = ^^^^^ Q"^(s, i, a, a; s, i, a, a) and iV/ = 2 
flavors. The coordinate n is subdivided into a spatial 
s and a temporal part t. By introducing energy eigen- 
functions one can show that these correlation functions 
decay exponentially in time, with the particle mass being 
the decay constant C{At) ^ exp(— mAt). On a toroidal 
lattice however with temporal extent T this exponen- 
tial decay appears as a cosh, C{At) ^ exp(— mAi) -|- 
exp(m(Ai — T)). Local masses m can be retrieved for 
every value of At by solving the implicit equations 

C{At -f 1) exp(-TO(A< + 1)) + exp(TO(Ai + 1 - T)) 



C{At) 



exp(-?7iAt) + exp(TO(At - T)) 



(28) 



with respect to m. Plateaus in the time dependence of 
the local masses m{At) exhibit the masses of the parti- 
cles. 



B. Smearing 

The extension of the local mass plateaus can be in- 
creased by enhancing the overlap with the ground state 
correlation function. We follow Ref. [|3| and apply the 
smearing matrix S = {F)^ , with 
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FIG. 17. Set of two-loop functions T'{At) on smeared 
sources and sinks plotted versus the cutoff 1. As in Fig. 
At increases as one steps down from the top curve. 



F(n, a, a; m, (3, b) 



1 



1 + 



S{n,a,a;m,P,b) (29) 



+1 ^ (yn,^,ia; b)S{n + fi; 
+Ul_f,^f,ia;b)S{n - fi;m)j 



m) 



choosing fc = 50 and I = 4. 

Source and sink smearing are readily accomplished by 
the replacements 

V', ^ r^ = (30) 

in the spectral propagator representation 



m{n,a,a;m,f3,b) = 

1 \^pl{n,a,a)}{iljl{m,/3,b)\ 



E 



A,: 



(31) 



Throughout this section we will assume the eigenmodes 
to be enumerated according to Eq. (pi]). 



C. Pion correlator 

The spectral representation of C7r(At) reads 

1 it)\^Pj {t)){i^j {t + At) I?/;, {t + At)) 



C^At) = 



XiXj 



h3,t 



{V'ilV'i){^j|V'i) 



(32) 



2J 



where we suppressed the brackets that represents the av- 
erage over the gauge fields. 

For C,r no early saturation over the entire At-range 
can be expected, as one can see by integrating Eq. (32) 
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FIG. 18. Two-loop correlator (Eq. estimated with TEA and SET, for the lightest quark mass and local sources and 
sinks. 



5^C.(Ai)=^-2 . (33) 

At I * 

Since all the contributions on the r.h.s. of Eq. (^) are 
positive the series is monotonically increasing. There- 
fore - contrary to the case of the topological charge - 
neither cancellation effects nor early saturation can be 
expected in this global quantity. But what about the 
regime of infrared physics described by the correlator, 
i.e. its asymptotic behavior in i? 

Let us consider the truncated spectral correlator 

I 

Ct(At)= ^ vi;(^,j- Ai). (34) 

In order to demonstrate the low-lying eigenmode domi- 
nance at large time separations we present in Fig. ( p^ ) 
a family of curves, C4(At), for the various time slices 
At, plotted against the spectral cutoff, at the light- 
est SESAM quark mass. It is gratifying to find that 
C^(At) for At > 7 shows a flat behavior in the regime 
I > 100. On the other hand for small time separations 



higher eigenmodes continue to add - in accordance with 
the idea of excited state contaminations. 

It is interesting to carry out a direct comparison of 
TEA with the standard propagator as computed by linear 
solvers on a local source, in order to see saturation occur 
in the region of interest, see Fig. (p^). We find good 
agreement in the asymptotic regime, 7 < At < 25. We 
also show the fit (fit range: [8, 15]) to the data from the 
inverter to the usual cosh parameterization: 

C^{At) = A cosh [m^{At - T/2)] . (35) 

A much more sensitive test of TEA is to look at local 
effective masses. In Fig. (|l^) we compare, at the cutoff 
value / — 300, the TEA results with the ones from stan- 
dard propagator analysis, for the lightest SESAM quark 
mass. The effects of smearing and varying quark masses 
are displayed in Fig. (|l5|), again for I = 300. As antici- 
pated, we do observe a clear tendency for improvement in 
the spectral approach with decreasing quark mass. Yet 
there remains a marked oscillatory behavior over the en- 
tire SESAM range of quark masses. Moreover we notice 
that smearing slightly improves the signal. 
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FIG. 19. Local 77' masses from TEA on local (LL) and smeared (SS) sources and sinks with ground state projection of its 
connected piece (OLGA). For comparison the n mass as obtained from the standard method is plotted. 



A synopsis on the cutoff dependence of ■m{At) is pre- 
sented in Fig. (^6|), as obtained on a particular configura- 
tion at the lightest quark mass. To avoid cluttering of the 
data we have spread out the different curves by means of 
the variable I At = 600 Ai -I- /. This survey plot is meant 
to convey an idea how the oscillation will dampen out 
with increasing cutoff I. 

Thus we have demonstrated that in the sea quark mass 
regime of the SESAM configurations, there is insufficient 
dominance of the low-lying eigenmodes to utilize TEA 
for a sensible calculation of the vr mass. 



D. ry'-correlator 

Let us consider next the flavor singlet pseudoscalar 
channel with the ground state particle 77'. As described 
in Eq. (|2^) C^/ differs from by the two-loop correlator 
T: 



C^,{At)^C\{At)-NfT{At) , 
the spectral representation of which reads 



(36) 



T{At) = 



(37) 

1 {tl:.j{t + At)\iPj{t + At)) 



A, 



where we suppressed again the brackets indicating the 
average over the gauge fields. 

Again we sum T[At) over At in order to learn about 
tea's potential in the two-loop situation 



(38) 



This is just the square of the 'topological charge', see 
Eq. (|l|). Therefore we might expect TEA to work as 



well as in section III 



We define the truncated two-loop correlator T' through 
^.A.^_ (39) 



T\At) = 



I 



1 {^,{t)\^,{t)) ^ 1 {^j{t + At)\iP,{t + At)) 



11 



0.55 

0.5 

o 0.45 
w 

ns 

^ 0.4 

ro 
o 

° 0.35 
0.3 
0.25 



SETk=0.156SS — X- 
SETk=0.1575SS □ 
TEASS - Q- 
71 mass ■ * 

e >© » % 



1=300 



X3 



X 



Local Masses 



6 8 

At 



10 



12 



FIG. 20. Comparison of the local OLGA 77' masses from 
TEA and from Stochastic Estimations on smeared sources 
and sinks. 



and plot in Fig. (17) the dependence of T'(At) on the 
cutoff I. Contrary to the pion propagator (see Fig. (p^), 
Fig. (|l^ indeed reveals good saturation of the spectral 
representation by the low-lying eigenmodes over the en- 
tire Ai-range. 

As a check for consistency we compare the local two- 
loop correlators from TEA and standard SET at our 
lightest quark mass in Fig. (|l^). They are seen to agree 
very well within their errors. Note that the TEA data 
show a much smoother behavior in At. Additional smear- 
ing for SET diminishes those fluctuations. We notice that 
TEA and SET data bear errors of equal size. We view 
this as an independent confirmation of previous claims 
that the errors on the 77' mass from state of the art SET 
analyses are dominated by gauge field noise |^J^. 

Next we come to the more stringent test: local effec- 
tive 77'-masses. According to Eq. (p6|), the ry'-propagator 
is the difference of one- and two-loop corrrelators, 
and T. The ground state contribution to the former, 
C^, can be determined very accurately by the standard 
methods (by iterative solvers) known from the octet spec- 
trum [ p2[ . Hence, it appears very natural to replace 
the one-loop correlator by its ground state component, 
C^(At), see Eq. (p5|). In the following, we perform a 
'one-loop groundstate analysis' (OLGA) by the extrac- 
tion of local masses, TO^/(Ai), from the combination 



C,,(Ai) = CS(At)-iV/T(At) 



(40) 



The results are presented in Fig. (|l9|), both for the 
lightest and heaviest sea quark masses of SESAM, with 
and without smearing. We find striking plateau forma- 
tion from the very first time slice onwards. 

As a consistency check and first test of the synthetic 
data approach, Eq. (|40|), we compare the local effective 
masses from TEA and SET in Fig. (|2^). The data points 
are seen to agree very well with each other, the TEA 
points being slightly less fluctuating. The horizontal lines 
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FIG. 21. Dependence of the local -q' masses m(At) for 
Ai = 1, 2, . . . 5 on the spectral cutoff. 

in Fig. (|l^) and Fig. (|2^) refer to the fitted plateau values 
for the 7]' masses. 

As ye t another test on systematic errors we plot in 
Fig. (|2l| ) - again for local and smeared wave functions 
and the lightest and heaviest quark masses - the depen- 
dence of the local masses m,,' (At) from TEA on the spec- 
tral cutoff I. It appears that the systematic errors from 
this cutoff are well under control, once we truncate the 
spectral representation with I ~ 150 and higher. Fur- 
thermore the data appear to support the idea that TEA 
improves when decreasing the quark mass. 

Finally we address the question, to what extent the 
r]' mass is influenced by the topological content of the 
configurations. By applying the cuts in Q, we subdivide 
the gauge field ensemble for the two K-values into two 
subsets each, with Q determined as described in section 
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FIG. 22. Dependence of the local OLGA r/' mass on the 
topological charge work similarly well for the lightest and 
heaviest sea quark mass of SESAM. 



III. The cuts are chosen such that each subset consists 
of about 100 configurations. The results from OLGA can 
be seen in Fig. (|2^). They clearly confirm the previous 
finding |^,^ , that topologically nontrivial gauge config- 
urations are the origin for the large 77' mass |34|,|35|] . 



V. DISCUSSION AND SUMMARY 

We have presented and validated a method to com- 
pute fermion loops from the low-lying eigenmodes of the 
hermitian form, Q = 75M, of the standard Dirac- Wilson 
matrix, M, in accordance with the expectation that these 
modes contain the essential physics associated with topo- 
logical fluctuations. Our truncated spectral approach 
(TEA) to is viable in the sense that it renders satis- 
fying results in the quark mass regime of state-of-the art 
fuU QCD simulations like SESAM on the basis of O(IOO) 
modes only. TEA has been verified both with respect to 
the topological charge and the two-loop correlator enter- 
ing the rj' propagator. 

The early onset of saturation for Tr Q^^ could be at- 
tained by proper ordering and subsequent partial sum- 
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FIG. 23. Pion correlation function from TEA with 600 
eigenmodes of AI on one configuration for k = 0.1575. 
One observes strong deviations from the standard propaga- 
tor (from linear solvers) determined on the entire gauge field 
ensemble. 



mations, configurationwise adapted to achieve cancella- 
tions from positive and negative eigenvalues. In this way 
the bulk of the higher mode contributions were shown to 
vanish. For the case of the 77' an early plateau formation 
of the local masses could be obtained by ground state 
projecting the connected piece of its propagator. 

We confirm previous results in the intermediate quark 
mass regime that the actual bottleneck in increasing ac- 
curacy of loop estimates is given by the gauge field noise, 
i.e. by the present limitation in ensemble sizes of QCD 
vacuum configurations. 

We found that the amount of work, i.e. the number of 
matrix vector multiplies, TV™""™^ needed per configura- 
tion in TEA compared to SET ^ for the lightest SESAM 
quark mass is roughly in the same ballpark^. This 
looks promising for the upcoming era of Teracomputing, 
where we shall deal with larger lattices and smaller quark 
masses. The reason is, that the Arnoldi method does not 
lose efficiency when entering deeper into the critical chiral 
regime, in contrast to Krylov solvers used within stochas- 
tic estimator algorithms which will suffer in convergence 
rate. 

Let us finally comment on the viability of the spectral 
approach applied to M instead of Q. To this end we com- 
pare in Fig. (|2^) the pion correlation function as deter- 
mined from 600 low- lying eigenmodes of M, see Fig. (|24|), 
with the one obtained from the standard method (solv- 
ing linear systems). Obviously, for the case of k = 0.1575 
which corresponds to a pion mass of the order of 730 
MeV, the chosen low-lying modes of M contain less in- 



"Actually the ratio iVrBA/^^SBT - 1-5 when using 300 
stochastic source vectors and 300 eigenmodes. 
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FIG. 24. The 600 calculated eigenvalues of M that enter 
Fig. ra. The curve represents a circle around the origin. 



formation about this function than the eigenmodes of Q. 
In order to trace down this discrepancy between M and 
Q we also carried out a full diagonalization on a quenched 
4^ lattice at /3 = 5.0 for the two cases. In Eig. (|2^), we 
plotted the resulting C.n-(At) versus the cutoff. While Q 
again shows a quite stable behavior similar to Eig. (12), 
the spectral approximation for M yields a very ragged 
cutoff dependence that requires all 3072 eigenmodes for 
a satisfactory representation of the correlator. We con- 
clude that the eigenmodes of M, being non orthogonal, 
suffer interferences among each other. Thus, for the 4^^ 
test case and for the sea quark masses used in the SESAM 
configurations, we have not been able to identify a lim- 
ited number of dominating eigenmodes. 

A detailed analysis of the 77' mass, based on OLGA and 
additional SET data, will be presented in a forthcoming 
paper ||6|. 
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